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Viscosity and thermal conductivity effects at first-order phase 
transitions in heavy-ion collisions 
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Kashirskoe sh. 31, Moscow 115409, Russia 
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Effects of viscosity and thermal conductivity on the dynamics of first-order phase 
transitions are studied. The nuclear gas-liquid and hadron-quark transitions in 
heavy- ion collisions are considered. We demonstrate that at non-zero thermal con- 
ductivity, K 7^ 0, onset of spinodal instabilities occurs on an isothermal spinodal 
line, whereas for k = istabilities take place at lower temperatures, on an adiabatic 
spinodal. 



I. INTRODUCTION 

There are many phenomena, where first-order transitions occur between phases of differ- 
ent densities. In nuclear physics various first-order phase transitions may take place in the 
Early Universe, in heavy-ion collisions and in neutron stars (e.g., pion condensation, kaon 
condensation, deconfinement and chiral phase transitions), see reviews IllQ. At low col- 
lision energies the nuclear gas-liquid (NGL) first-order phase transition is possible 

m. 

At high collision energies the hadron-quark gluon plasma (HQGP) first-order transition may 
occur (see e.g. Ref. [2|). Within a hydrodynamical approach dynamical aspects of the NGL 



7 



llj . An important role of effects of 



and HQGP transitions were recently studied in Refs. 

non-ideal hydrodynamics was emphasized. In this talk, we review some results obtained in 
Refs. 
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II. GENERAL SETUP 

We are interested in a description of long wavelength phenomena at a first-order phase 
transition. Thereby we accept an assumption that the velocity of a fluctuation (seed) u is 
much slower than the mean thermal velocity. Then a description is possible in the framework 
of the standard system of equations of non-relativistic non-ideal hydrodynamics: the Navier- 
Stokes equation, the continuity equation, and equation for the heat transport: 
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mn[dtUi + {uV)ui] = -V^P + 



V (^kUi + ^iUk - -^Sikdivuj + C^ifcdivu 



1) 



T 



dtn + div{nu) = 0, (2) 

ds ^ ( 2 \^ 

— + div(su) = div(KVT) + 77 \Vi,Ui + V - ^^ifcdivuj + C(divu)^ , (3) 



n is the density of the conserving charge (here the baryon density), m is the (baryon) mass, 
P is the pressure, 77 and C are the shear and bulk viscosities, d is the dimensionality of space, 
T is the temperature, s is the entropy density, k is the thermal conductivity. 

The simplest (yet non-trivial) example illustrating principal features of a first-order phase 
transition in a mean-field approximation is the Van der Waals fiuid. The adiabatic trajec- 
tories, s = s/n ~ const, for an expansion of a uniform fireball to vacuum, are shown in 
Fig. 1. The super-cooled vapor (SV) and the overheated hquid (OL) regions are between 
the Maxwell construction (MC) and the isothermal spinodal (ITS) curves, on the left and 
on the right respectively. The adiabatic spinodal (AS) curve bounds the AS region from 
above. For Scr > s > smc2, where corresponds to the value of the specific entropy s at 
the critical point and the line with Smc2 passes through the point n/ucr = 3 at T = 0, the 
system traverses the OL state (the region OL in Fig. 1), the ITS region (below the ITS line) 
and the AS region (below the AS line). For s > Scr the system trajectory passes through 
the SV state (the region SV in Fig. 1) and the ITS region. 

All thermodynamic quantities can be expanded near an arbitrary reference point (n^, T^), 
or (nr,Sr). Considering the problem in {n,T) variables it is convenient to take the reference 
point in the vicinity of the critical point {ncr,Tcr) but outside the fluctuation (critical) 
region {ricr T Sn^, Tcr ^ 5T^) assuming that the latter is very narrow. Note that even in the 
fluctuation region, the mean-fleld treatment can be used provided one considers the system 
at time scales, being shorter than the scale responsible for a development of long-scale critical 
fluctuations. Thus we further put (rir,?!) = {ricr-iTcr)- 



3 



Let us construct a generating functional in the variables Sn = n — Ucr, ST = T — T^r, also 
called the Landau free energy, such that 6{6FL)/6{6n) = P — Pj + Pmc- 

where e = Pf — Pmc is expressed through the pressure at the MC. The maximum of the 
quantity e is e"* = 4Af ^/(3-\/3)- For the sake of convenience, we will use e normalized to its 
maximal value, 7^ = le/e*"! where < 7^ < 1. The first term in Eq. (jl]) is due to the surface 
tension, SFL^snrf = ^''S', 5* is the surface of the seed. For the Van der Waals equation of state: 

Using Eq. ([2]), we rewrite Eq. ([T]) in the dimensionless variables 6p = vip, = Xi/l, 
i = 1, ■ ■ ■ , d, T = t /to as 

= A, (a,^ + 2^(1 -^')+^-^)^ (6) 



(i = 2{d — l)/d. Only linear terms in the velocity u were kept in deriving this equation. 
Since f ^ oc —6T, processes in the vicinity of the critical point are proven to be very slow. 
This is known as the critical slowing down phenomenon. 

Eq. ([6]) should be supplemented by Eq. (|3]) for the heat transport, which owing to Eq. ([2]) 
after linearization simplifies to 



dtSs - s,{n,)-^dt6n] = k^AST. (7) 



The variation of the temperature is related to the variation of the entropy density s[n, T] by 
6T ~ Ti.(cy,r)~^ i^s — {ds /dn)T^r5n) , where cy,r is the specific heat density. 

Note that Eq. ([6]) differs from the standard Ginzburg-Landau equation exploited in phe- 
nomenological approaches. The difference disappears, if one sets the bracketed term in the 
r.h.s. of Eq. (jG]) to zero. Superficially this simplification is legitimate, if space-time gradi- 
ents are small. However, for a seed prepared in a fluctuation at t = with a distribution 
5p{t = 0,r), the initial condition ^^^qI'^ \ t=o — should be fulfilled. Otherwise, because of 
a positive kinetic energy contribution the probability of fluctuations is suppressed. On the 
other hand, the above initial conditions cannot be simultaneously fulfilled, if a differential 
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equation is of the first-order witli respect to time derivative. Therefore, there exists an 
initial stage of the dynamics of seeds (t < tinit), which is not described by the standard 
Ginzburg-Landau equation. 

The time scale for the relaxation of the density following Eq. (jG]) is tp oc R, where R is 
the size of a seed (as we will show below, at some stage overcritical seeds grow with constant 
velocity), and the time scale for the relaxation of the entropy/temperature, following ([7]), is 

tT = R^Cv,r/Hr OC R^ . (8) 

The evolution of a seed is governed by the slowest mode. Thus, for tT{R) < tp{R), i.e for 
R < -Rfog (-Rfog is the typical seed size at which tp = t^), dynamics of seeds is controlled by 
Eq. ([6]) for the density mode. For seeds with sizes R > -Rfog, oc R? exceeds tp oc R and 
growth of seeds is slowed down. Thereby, the number of seeds with the size R ~ Rfoo- may 

un 

increase with time. Estimates p, 18| show that for the HQGP phase transition Rfog ~ 0.1 — 1 
fm and for the NGL transition i?fog ~ 1 — 10 fm < Rfh{tf,o.), where Rfb{tf,o.) is the fireball size 
at the freeze out, t/ o. is the fireball evolution time till freeze out. Thus, thermal conductivity 
effects may manifest themselves in heavy-ion collision dynamics. 

Note that seeds of a new phase are produced in an old phase owing to short-scale fluc- 
tuations. Fluctuations are not incorporated in the above hydrodynamical equations defined 
in terms of mean-field variables. Contributions of short-scale fluctuations can be simulated 



by a random force induced in Eqs. (E]), ([7]), cf. |12|. 

There are only two dimensionless parameters in Eq. ([6]), e and /3. The parameter e is 
responsible for a difference between the Landau free energies of the metastable and stable 
states. For tp ^ tx (isothermal stage), e ^ const and dependence on the latter quantity 
disappears because of A^e ^ 0. Therefore, dynamics is controlled only by the parameter 
P, which characterizes inertia (enters together with the second derivative in time). This 
parameter is expressed in terms of the surface tension and the viscosity as 

/3 = i32T,r)-'[dr], + CVa^o^- (9) 

The larger viscosity and the smaller surface tension, the effectively more viscous (inertial) 
is the fluidity of seeds. For /3 ^ 1 one deals with the regime of effectively viscous fluid 
and at (3^1, with the regime of perfect fluid. For the NGL phase transition we estimate 
(3 ~ 0.01. For the HQGP phase transition /3 ~ 0.02 — 0.2, even for the very low value of 



the ratio rj/s ^ l/(4vr) {3, Q]- Thus one deals with effectively very viscous fluidity of density 
fluctuations in both NGL and HQGP transitions. 

III. DYNAMICS OF SEEDS IN METASTABLE AREA 

Let us consider the stage tp 3> tr- For an expanding system we have to assume that 
typical time for the formation and evolution of a fluctuation of our interest t/orm + tp is 
much smaller than the typical flreball expansion time t/.o.. Let us consider the situation, 
when at very slow expansion with s{t) ~ const spatially quasi-uniform spherical flreball of a 
large radius Rfhif) enters either the OL- state or the SV- state (see the corresponding curves 
in Fig. 1). In case \n—ncr\/ncr ^ 1, i.e. in the vicinity of the critical point (rzcr, Tc^), solution 
(|6]) describing dynamics of the density in the fluctuation is presented (in the dimensional 
form) as [3, Q : 



on[t, r) ~' 



m 



±th- 



{6n)cor, (10) 



where the upper sign corresponds to the evolution of bubbles and the lower one to the 
evolution of droplets, d = 3, and the solution is valid for |e/ {XcrV^{T)) \ ^ 1. The correction 
{6n)cor is responsible for the exact baryon number conservation. Considering r in the vicinity 
of a bubble/droplet boundary we get equation describing evolution of the seed size: 

Pt^d'^Rn _ 3e 21 todRn , . 

~Y df^ ~ 2XcrV^{T) ~R^~ l~dr' ^ ' 

Following this equation a bubble of an overcritical size R > R^r = 'ilXcrv'^ {T) / {"it) of the 
stable gas phase, or respectively a droplet of the liquid phase, been initially prepared in a 
fluctuation, will grow. On the early stage of the evolution the size of the bubble/droplet 
Rn{t) > Rcr grows with an acceleration. Then it reaches a steady grow regime with a 
constant velocity Uas = a v»'^lT)to °^ le\Tcr — T\^/'^. In the interior of the seed 6n ~ ^v{T)/m. 
The correction {6n)cor — 'v(T)Rl^(t) / {mR^j-j^) is very small for Rnit) ^ Rfbif). 
Substituting Eq. ffTOj) to Eq. ([7]) for T = const we obtain 

r-Rn{t) , e " 



±th- 



+ {6n),or > . (12) 



/ 2XcrV%T)_ 

While the temperature is constant in the interior and exterior of the seed, the entropy and 
the density are different in the interior and exterior regions. Nevertheless the approximation 
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of a quasi-adiabatic expansion of the system might be used even, when the system reaches 
metastable region, provided the gas of seeds is rare, or v{T) is small. 

In Fig. [2] we demonstrate numerical solutions of the hydrodynamical equations in two di- 
mensions {d = 2) at the stage tp ^ tr- We take T/Tcr = 0.85 and compute the configuration 
for 7] ~ 45MeV/fm^ and for (3 = 0.2 (effectively large viscosity). The choice Tcr = 162 MeV, 
"n/nsat = 1-3 {usat = 0.16 fm~^) is relevant for the HQGP phase transition. We see that un- 
dercritical seeds (disks) dissolve rather rapidly (typical time is ~ several fm) but overcritical 
seeds grow very slowly. Similar solutions exist for bubbles. Therefore, one can hardly expect 
to observe a manifestation of large size droplet/bubble remnants in heavy- ion collisions. 

The limit k = is specific. In case \n — np^maxl/np^max 1, i-e. in the vicinity of the 
point {np^max, sp^max), Corresponding solutions of Eq. (I6j) describing dynamics of the density 
can be presented in the form ffTOj) with the only difference that 6T should be replaced by 6s, 
\cr by Xp^max and v(T) by v{s). Dynamics of -R„(t) is determined by Eq. (fTTj) . where one 
should replace values calculated at {ricr^Tcr) to the corresponding values at (np,max, sp^maxO- 
From Eq. ([6]) we obtain 55 = = {dsnp^max — s p^max^n) / n^p with 5n given by Eq. (|TOl) . 
From known values 5s and 5n we can define 5T . Therefore, not only the density, n, but 
also the entropy density, s, and the temperature, T, vary in the surface layer and exhibit 
different values inside and outside a seed. Contrary, the value s remains constant. 



IV. INSTABILITIES IN SPINODAL REGION 



In this section the "r" -reference point can be taken arbitrary, therefore, we suppress the 
subscript "r". To find solutions of the hydrodynamical equations we put, cf. [3, 0], 

5n = 6nQexp['jt + ^pr], 6s = 5soexp[7t + ipr], T = + 5Toexp[7t -|- ipr], (13) 

where T> is the temperature of the uniform matter. From the linearized equations of non- 
ideal hydrodynamics we find the increment, 7(p), 
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-P 



mn 



(14) 



1 + Kp'^/{cvl)_ 

u\- = m~^{dP/dn)'p and u| = m~^{dP/ dn)^ are speeds of sound at constant temperature 
and entropy, respectively. Eq. ( !T^ has three solutions. Expanding the solutions for small 
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momenta (long-wave limit) we find 

7i,2 = ±iUsP + 



p2 



2 2 

^3 = 3— 



2 2 I 2 2 2 I ^' 

ujuj< \ ujcy mncvuj I 



(15) 
(16) 



MfCy 

The solutions 71^2 correspond to the sound mode in the long wavelength limit, whereas 73 
describes the thermal transport mode. Below the ITS line (and above the AS line) u% < 0, 
ii| > 0, solutions 7i^2 correspond to an oscillation and damping, whereas 73 describes an 
unstable growing mode. Below the AS line, since there < and < 0, the modes 
exchange their roles: the sound modes become unstable, while the thermal mode is damped. 

For sufficiently high thermal conductivity and not as small p, /tj9^/(cy I7I) ^ z/, z/ = 
(m| — «!.)/(— My), within the ITS region the most rapidly growing mode (7 ~ 7^, and 
p ~ pm) is the density mode: 

i2VP + l)idrj + Cy (2^/?+l)c- 

Using these expressions we may rewrite the condition of a high thermal conductivity as 
K. ^ ucyy/c- Typical radius of structures is Rm ~ l/Pmj it decreases with increasing \dT\. 
The amplitudes of the temperature and density are related by 

. ,„„^*_^£^A^. (17) 

Cyn [1 + K/[^yCCv)\ 

Therefore, the assumption of spatial homogeneity of the system fails right after the ITS 
region is reached. An aerosol (mist) of bubbles and droplets is formed for a typical time 
taer ~ ^/im- the Van der Waals equation of state we find that STq/Suq > 0, i.e. 
the temperature is larger in denser regions. However, the amplitude of the temperature 
modulation is rather small for k ^ ucvy/c- 

Let us consider the case of zero shear and bulk viscosities and non-zero, but small thermal 
conductivity. For —u^ <^ 1, i.e. shghtly below the ITS line, we get 

^- = ^-"4^' Pl^-^T/i2c). (18) 

Therefore, in both considered cases the instability occurs at the ITS hne. However, for 
K ^ vcvy/c the most rapidly growing mode corresponds to 7^'^^ and in the opposite limit 
K < vcv^Jc, to 7^). 
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In Fig. |3]we show the time evolution of the density wave amphtudes given by the first 
Eq. f|T3l) . for an undercritical value of p, p < Pcr = V^/l (left panel), and for an overcritical 
value (right panel) for the same choice of the parameters T/Tcr = 0.85, Tcr = 162 MeV, 
fi/usat = 1-3, as in Fig. |2l We see that the evolution is more rapid compared to the 
configuration presented in Fig. |2]and the characteristic time scale tch ~ 10 fm is comparable 
with the fireball expansion time tf ^.- Thus we may conclude that in heavy- ion collisions 
during expansion of the fireball the system may linger in QGP phase at T < T^r- This 
means that the equilibrium value of the critical temperature of the phase transition might 
be significantly higher than the value which may be manifested in growth of fiuctuations in 
experiments with heavy ions. Fluctuations grow more rapidly with decrease of T below Tcr, 
tch ~ l-5Tcr./|(5T| fm. Far from the critical point rapid grows of fiuctuations reminds effect 
of a warm champagne. There are prospects for observing specific signatures of fiuctuations 
with a typical size, defined by 1/pm, in heavy-ion collisions. Fluctuations of this kind might 
be distinguishable from ordinary statistical fiuctuations. 

The limit k = is again specific. Below the ITS line and above the AS line, the thermal 
mode 73, which drives the system towards equilibrium for a small thermal conductivity, does 
not exist for k, = 0. Therefore, the evolution in the spinodal region is entirely governed by 
adiabatic sound excitations. The increment 7 is given by Eq. (HM with ut replaced by uf. 



2 2 
1 =-P 



2 , idr] + Ch , 2 

U z H h cp 

mn 



(19) 



Therefore, contrary to the case k 7^ 0, instability appears, when the system trajectory crosses 
the AS line rather than the ITS line. The value m| = —Xp^maxv'^{s) = 155sTpmaa;/(512m) 
for the particular case of the Van der Waals fiuid in the vicinity of the point {np^max, Tp^max) 
(taken as the reference point). This result holds also in the case of ideal hydrodynamics, 
where in addition to k = 0, the viscosity coefficients [r], () are zero. From Eq. ([7]), we find 

5s, = 6n,s/n, 5To = 5no^^i^^^^^^^^^^. (20) 

Thereby, the temperature is modulated similar to the entropy density. For the Van der Waals 
equation of state we find that STo/Sno > 0. The amplitude of the temperature modulation 
is larger than in case of /t 7^ 0, see Eq. fll7p . 

Concluding, for any k 7^ the solutions of Eq. ( !T^ result in the onset of the instability 
already for u"^ < (i.e. below the ITS line rather than below the AS line). Since in reality k 
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is indeed nonzero, spinodal instabilities start to develop when the trajectory crosses the ITS 
line rather than the AS, i.e. at significantly higher temperatures. This favors an observation 
of signals of the spinodal decomposition in the HQGP phase transition in heavy-ion collisions. 

We expect that owing to a manifestation of non-trivial fiuctuation effects (especially, 
of the spinodal decomposition at first-order hadron-quark transition) a non-monotonous 
behavior of different observables as function of coUisional energy may be observed. This 
experimental analysis will be possible at RHIC, FAIR and NICA. Owing to properties of 
spinodal decomposition a manifestation of specific structures with typical spatial size may 
serve as a promising signal of the QCD first-order phase transition. 
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FIG. 1: The phase diagram of the Van der Waals equation of state, T(n)-plane. The bold solid, 
dashed and dash-dotted curves demonstrate the boundaries of the Maxwell construction, the spin- 
odal region at T =const and s =const, respectively. The short dashed lines show adiabatic trajec- 
tories of the system evolution: the curve labeled Scr passes through the critical point; Sm, through 
the maximum pressure point P{np^max) on the P{n) plane. 




11 

2 
1.5 

u 

-5- 1 

Q. 

0.5 

% 0.1 0.2 0.5 1 1.5 0.5 1 1.5 
P/P r/L r/L 

cr 

FIG. 2: The isotherm for the pressure as a function of the density, with initial and final configu- 
rations shown by the dots (left column). The dashed vertical line shows the MC. The initial state 
represents the stable liquid phase disk (d = 2) in metastable SV. Middle column demonstrates the 
time evolution of the density of the overcritical liquid disk. The numbers near the curves (in L) 
denote time moments; r = + y^, \ST/Tcr\ = 0.15, L = 5 fm. Right column, the same for the 
undercritical liquid disk. 




12 




FIG. 3: Time evolution of the wave amplitudes f{t) defined as 6n{t) normalized to the amplitude 
of the initial disturbance. Solid line is for effectively small viscosity (/3 = 10) and dash line, for the 
large viscosity (/3 = 0.1). Left panel: the undercritical wave number p = 2/L (growing modes). 
Right panel: the overcritical value p = 8/L (oscillation modes for large /3 and damped modes for 
small /3). Other parameters are taken to be the same, as in Fig. [2j 
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FIGURE CAPTIONS 

1. The phase diagram of the Van der Waals equation of state, T(n)-plane. The bold solid, 
dashed and dash-dotted curves demonstrate the boundaries of the Maxwell construc- 
tion, the spinodal region at T =const and s =const, respectively. The short dashed 
lines show adiabatic trajectories of the system evolution: the curve labeled Scr passes 
through the critical point; s^, through the maximum pressure point Pijip^max) on the 
P{n) plane. 

2. The isotherm for the pressure as a function of the density, with initial and final con- 
figurations shown by the dots (left column). The dashed vertical line shows the MC. 
The initial state represents the stable liquid phase disk {d = 2) in metastable SV. 
Middle column demonstrates the time evolution of the density of the overcritical liq- 
uid disk. The numbers near the curves (in L) denote time moments; r = 
\6T/Tcr\ = 0.15, L = 5 fm. Right column, the same for the undercritical liquid disk. 

3. Time evolution of the wave amplitudes f{t) defined as Sn{t) normalized to the ampli- 
tude of the initial disturbance. Solid line is for effectively small viscosity (/3 = 10) and 
dash line, for the large viscosity (/3 = 0.1). Left panel: the undercritical wave number 
p = 2/L (growing modes). Right panel: the overcritical value p = 8/L (oscillation 
modes for large /3 and damped modes for small /3). Other parameters are taken to be 
the same, as in Fig. [2l 



